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Summary 


Until  recently,  computations  of  space  plasma  flow  over  a  spacecraft  have  been 
unstable  for  ratios  of  spacecraft  dimension  to  Debye  length  typical  of  the  low  earth 
orbit  environment.  We  present  calculations  of  the  spacecraft/environment  inter¬ 
action  based  on  two  computer  codes,  MACH  and  POLAR. 

We  have  developed  MACH,  an  inside-out  particle  tracking  code,  for  the  purpose 
of  validating  the  physics  of  POLAR  in  regimes  where  there  are  no  comprehensive 
theoretical  or  experimental  results.  While  the  spacecraft  which  can  be  treated  by 
MACH  are  restricted  to  simple  geometries,  the  methodology  is  more  fundamental 
than  POLAR.  M.\CH  generates  self-consistent  solutions  within  the  context  of 
quasisteady  Viasov  plasma  flow  and  achieves  Debye  ratios  previously  unobtainable. 

POLAR  uses  a  three-dimensional  finite-element  representation  of  the  vehicle 
in  a  staggered  mesh.  The  plasma  sheath  is  modeled  by  outside-in  particle  tracking. 
Solutions  for  the  plasma  flow,  wake  and  vehicle  charging  are  obtained  by  Vlasov- 
Poisson  iteration;  charge  stabilization  techniques  make  the  results  virtually  insensi¬ 
tive  to  the  Debye  ratio.  P01..\R  reproduces  the  Laframboise  static  pl.asma  solutions 
for  spherical  probes  and  fits  the  Makita-Kuriki  probe  data  for  spheres  in  a  flowing 
plasma  in  regions  where  comparisons  are  valid. 

POLAR  and  ^L\CH  solutions  for  the  particle  and  electrostatic  potential  structure 
of  the  wake  of  a  charged  disk  in  a  low-altitude  flow  are  shown  for  Macli  numbers 
4,  a,  and  H.  We  compare  the  results  of  the  codes  to  each  other  and  to  simple 
analytic  approximations  as  part  of  an  effort  to  compute  the  nature  of  the  spacecraft 
wake  and  to  establish  the  validity  of  the  codes. 
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Computer  Models  of  the  Spacecraft  Wake 


1.  INTRODU(;TlO^ 

Knginee ring  design  of  spacecraft  requires  that  computer  codes,  such  as  I'Ol.AK, 
be  developed  to  simulate  the  spacecraft/environment  interaction,  in  order  to  test 
such  codes,  there  arc  three  general  lines  of  attack  available;  comparison  with 
tl)  experiment,  (2)  analytical  results,  and  (3'  other  computer  codes.  These 
approaches  are  limited  bv  the  fact  that  comprehensive  experimental  results  are  not 
yet  available  and  that  the  interaction  between  ;i  spacecraft  and  a  flowing  plasma  is 
sufficiently  complex  that  no  comprehensive  exact  solutions  exist  for  conditions 
typical  of  the  low  llarth  orbit  environment.  .As  a  consequence,  the  most  effective 
too!  fo'-  code  verification  is  comparison  with  other  computer  codes.  .As  part  of 
the  effort  involved  in  the  development  and  verification  of  the  1H)1.;\K  code.  We  liave 
developed  an  independent  compute  r  code,  called  AI.ACll,  designed  to  validate  the 
phvsical  and  numerical  methods  of  l'()l,.\H.  M.ACll  is  constructed  essentially  from 
fust  [irinciples,  and  can  be  used  to  examine  the  fundamental  plasma  processes 
in  tlu'  spac:eci’aft  inter.iction  and  the  formation  of  the  wake.  The  purpose  of  this 
repoi't  IS  to  disv  uss  the  results  ol  the  M.AC  ll  and  l'()l..Al{  computer  models  of  the 
sp.u  ei- raft  cn  vi  ronment  interaction  and  wake. 

I  lieceived  for  publieation  23  .lulv  If'Ttil 
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The  methodologies  of  POLAlf  and  MACH  are  somewhat  different.  It  is  these 
differences  that  make  possible  the  validation  of  one  code  by  another.  I'tfLAH 
(Potential  of  l.ai-ge  Spacecraft  in  the _Auroral  Hegions)  is  designed  as  an  engineer¬ 
ing  tool  to  compute  the  plasma  interaction  of  largo  electrically  complicated  space- 
ci'aft  in  low  Earth  orbit.  The  essential  differences  between  PtfLAH  and  its  prede¬ 
cessor,  NASC.AP  (which  was  designed  for  the  geosynchronous  environment),  are 
that  the  effects  of  short  Debye  lengths,  energetic  auroral  electrons,  and  the  geo¬ 
magnetic  field,  along  with  the  formation  of  a  supersonic  ion  wake,  are  taken  into 
account.  .-\n  essential  constraint  on  POLAK  is  that  the  physical  and  numerical 
algorithms  used  must  allow  computations  to  be  completed  in  a  reasonable  length  of 
time,  for  example,  on  tlie  order  of  hours  or  less  per  solution.  To  meet  this  con¬ 
dition,  a  number  of  novel,  but  essentially  untested,  physical  assumptions  have  been 
incorporated  into  the  code.  s  a  consequence,  the  results  are  to  be  regarded  as 
a  working  model  whose  validity  is  to  be  tested.  On  the  other  hand,  PCJl.AH  is 
capable  of  accounting  for  the  full  complexity  of  the  three-dimensional  interaction 
witliin  the  context  of  the  assumptions  made. 

The  methodology  of  MACH  (Mesothermal  Auroral  Charging)  is  somewhat  closer 
to  first  pr.nciples  than  POLAR.  MACH  solves  the  same  Vlasov-Poisson  equations 
ns  POLAR  but  under  different  conditions.  The  most  important  difference  is  that 
no  assumption  about  the  sheath  edge  is  made:  the  sheath  and  the  wake  are  deter¬ 
mined  as  the  self-consistent  solution  to  the  \'lasov-Poisson  equations.  However, 
the  geometry  of  M.ACH  is  restricted  to  two  dimensions,  precluding  detailed  appli¬ 
cation  to  real  spacecraft.  In  addition,  MACH  is  restricted  to  a  single  Maxwellian 
plasma  and  fixed  surface  potentials  (that  is,  the  surface  potentials  ai-e  given  and 
not  determined  as  a  solution  to  a  boundary  value  problem).  MACH  pays  a  steep 
price  for  its  adherence  to  fundamentals;  it  is  much  less  robust,  requires  a  great 
deal  of  operator  attention  and  requires  about  an  order  of  magnitude  more  execution 
time.  The  M.\CH  program  is  an  adaptation  of  TD\V.-\KE  (Parker'). 

In  describing  and  comparing  the  two  codes,  we  look  first  at  the  Poisson  solvers. 
M.ACH  uses  a  first  order  successive  over-relaxation  method,  while  PCfLA  R  uses 
.1  more  sophisticated  second  oi'der  finite  element  scheme  with  a  conjugate  gradient 
solver.  Both  codes  have  passed  a  variety  of  l.aplacian  (zero  plasma  density)  and 
nonflowing  Poisson  (nonzero  plasma  density'  tests.  These  tests  include  the 
I  ..if  l  amboise  static  plasma  solutions  for  spherical  probes  and  tlie  Makita-Kuriki 
probe  data  for  spheres  in  a  flowing  plasma  in  regions  where  comparisons  are  valid, 
rtius,  we  currently  di.snnss  Poisson  solver  diffe ri'nres  as  a  source  of  discrepanrv 
in  tlic.se  comparisons. 

1.  Parker,  L.  U  .  (IhTti)  Conti’actor  Report  .No.  CK-144159. 


At  this  stage,  we  narrow  our  study  to  three  possible  sour-ces  of  discrepancy; 

(1)  the  \’lasov  solution  for  the  number  density,  (2)  the  methods  of  stabilizing  tire 
Poisson  solvers,  and  (3)  grid  resolution  effects.  The  fii’st  is  the  greatest  physical 
interest,  the  second  is  of  numerical  origin,  while  the  third  is  not  tundamental  but 
must  be  dealt  with. 

The  two  codes  differ  the  most  in  their  Vlasov  solutions.  Tlie  inside-out  method 
used  in  \LAlT!  is  a  standard  mctliod  used  in  spacecraft  chai-ging  models.  .‘\t  eadi 
grid  point  in  the  problem,  each  trajectory  is  traced  backwards  in  time  to  discovi'i' 
whence  the  particle  came,  from  the  object  or  the  ambient  plasma.  Use  of 
Liouville' s  theorem  allow s  the  construction  of  the  distribution  function  at  the  point 
and  the  particle  number  densities  from  the  first  moment  of  the  dist labution  function. 
The  method  solves  the  \Tasov  equation  exactly  subject  only  to  numerical  appro.xi- 
mations,  but  demands  an  enormous  number  of  trajectories,  which  ai-e  subject  to 
cumulative  numerical  errors,  and  consequently  a  large  amount  of  computation  time. 
Tins  method  is  at  present  too  expensive  for  a  thi'ce-dimensional  code  placing  a 
priority  on  efficiency.  Thus  POI.AH  uses  an  approximate  method  that  is  efficient 
and  whose  accuracy  is  the  focal  point  of  this  comparison.  For  brevity,  let  us 
assume  a  convergent  sequence  of  \Tasov-Poisson  iterations.  Following  a  Poisson 
iteration,  the  solution  is  inspected  for  a  sheath  edge  (nominally  defined  by  electro¬ 
static  potential  -  0.-19  kT/e)  and  outside-in  particle  trucking  is  used  to  calculate 
the  attracted  specie  number  density  witliin  the  sheath  by  weighting  by  the  time  each 
particle  spends  in  each  volume  element.  Outside  the  sheath,  densities  of  attracted 
particles  are  calculated  from  the  geometric  shadow  of  the  object,  Hoth  W.VCll  and 
P(.)I..\1{  assume  Boltzmann  electrons,  valid  provided  that  the  electron  Mach  number 
is  small,  that  the  surface  potential  is  sufficiently  negative  that  electrons  do  not 
reacli  the  spacecraft,  and  that  no  potential  wells  develop.  Both  programs  take 
<I>  -  0  (where  <t>  is  the  electrostatic  potentiaB  as  the  boundary  condition  at  the  outer 
grid  boundary  and  start  from  an  initial  density  distribution  appropriate  for  tlie  flow 
of  a  collisionless  neutral  gas  past  the  spacecraft.  Successive  t  lasiw- Poisson 
iteration  cycles  are  computed  until  a  converged  solution  is  achieved. 

The  Vlasov-l'oisson  iteration  process  is  inherently  unstable  when  the  Debye 

a 

length  is  less  than  the  grid  spai-ing  i  Parker  .ind  Sullivan”).  To  date,  only  two  tech¬ 
niques  exist  for  treating  this  instability;  Pie.ird  iterativw  mixing  (Parker'  and 
.Sullivan  '  and  chai'ge  stabilization  (Cooke  et  al'  ).  The  latter,  which  is  used  in 

2.  Pai'ker,  I..\\.  .  .md  Sulliv.tn,  P.t’.  \.\S,\  liejioi't  .\o.  TN  -  D-TIO'', 

3,  Cooke,  D.  i,.  ,  Katz,  1.,  Manilell,  l.illi'V,  i.K.,  ,1 1'.  ,  .uid  Kubin, 

(19115)  .A  I'll  I'ee -dinrension.i  I  C.ilculation  of  Shutth-  t  harging  in  Polar-  (irbd. 
Spacecraft  Fnvi  i'onment:i  1  lute  I'actions  I'e,  hnology- 1  9ir(.  NASA  (  P-  235‘', 
Al-C'.l.-TK-Hri-OOlU,  p.  d'd". 
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POLAI{,  is  a  correction  appJied  at  each  iteration;  it  involves  a  charge  density 
adjustment  whenever  it  appears  to  be  numerically  amplified  by  finite  gi  id  resolu¬ 
tion.  The  iterative  mixing  method  pi-csumably  has  less  impact  on  the  physics,  but 
this  has  not  been  conclusively  demonstrated.  This  accuracy  issue  is  addressed  by 
using  both  methods  in  MACH. 

An  efficient  three-dimensional  code  must  keep  execution  time  within  reason¬ 
able  bounds.  In  order  to  attain  this  goal,  FCJI.AH  is  typically  run  with  minimum 
spatial  resolution,  so  that  the  effects  of  the  finite  gi'id  size  is  a  practical  concern 
and  must  be  dealt  with  in  order  to  develop  a  working  model. 

Our  ultimate  aim  is  to  develop  HOL.VH  to  calculate  the  spacecraft  interaction 
in  a  quite  general  way.  The  approach  which  we  will  present  here  is  a  bootstrap 
method;  begin  with  a  restricted  code  (M.ACH)  whicli  does  tlie  fundamental  physics 
correctly,  check  it  by  available  analytical  and  physical  means,  and  then  use  it  as 
the  standard  of  comparison  for  the  more  general  I'Oile  ( I’tJL.V H) .  In  this  report  we 
pi'esent  results  from  MACH,  possibly  the  first  coiie  to  i-esult  m  a  satisfactory  wake 
.md  inter.iction  model,  compare  them  with  analytical  results  wliere  applicable,  and 
discuss  the  POL.AU  code  in  comparison  with  MACH.  Our  conclusion  will  be  that 
M.'XC'H  and  POL.AU  represent  generally  accurate  qualitative  and  quantitive  models 
of  the  spacecraft  interaction  and  wake,  features  of  the  solution  wluch  aiv  not  yet 
understood  are  discussed. 

■A  word  about  the  nature  of  the  presentation  is  appropriate.  During  the  pi-epara- 
tion  of  tfie  report,  the  eomputer  system  which  runs  both  models  faijim'.  .As  a  re¬ 
sult,  the  discussion  is  based  on  two  sets  of  corriputations:  MACH  at  Mach  A  and  ii, 
on  the  one  hand,  and  POL, .AH  and  MACH  -it  M.icii  5  on  the  other,  the  ciiuice  of  Mach 
numbers  aepends  on  the  results  in  existence  at  the  time  of  the  failure.  U  hile  tlie 
report  could  have  been  based  on  Mach  a  alone,  the  Mach  number  dependence  was 
deemed  important  to  illustrate. 


THE  W  AKE  MODEL 


To  compute  models  of  the  spacecraft  interaction,  we  have  specialized  the  space¬ 
craft  in  both  MACHl  and  P(A[..AU  to  charged  disks  in  flowing  plasmas.  Siudi  a  prob¬ 
lem  in  an  unmagnetized  plasma  can  be  described  in  terms  of  three  dimensionless 
parameter.s  (assuming  that  the  relative  orientation  of  the  bulk  velocity  and  the  disk 
normal  is  fixed): 


-  e4.;  kT 


(!) 


=  a /A 


D 


(2) 


T'-y^.p  » l."  I  'J*  ■  IT  >  yj  ipj  w  w«.i  ir?  i 


L’nu^upma  w  v- 


M  =  \'  ;\  {  2  kT  '  m  '  .  ( ) 

1 

Here,  e  is  magnitude  of  the  electron  charge,  k  is  Holt/.maiin' s  constant,  T  is  the 

temperatuie  of  the  ambient  plasnui  (assumed  to  be  the  same  for  L'lecti-ons  and  ions), 

a  is  the  disc  radius,  'here  defined  us  the  ratio  of  thermal  speed  to  [ilasma  Ire- 

quencsi  IS  the  ambient  ion  Ifebye  U-ngtli,  M  -  is  tlie  ambient  Macli  number,  t  is 

the  ambient  bulk  \elocitv,  and  ni .  is  the  ion  niass;  tlie  ions  are  assumed  to  be 

1 

singly  loni/ed.  In  a  magnetised  plasma,  the  ion  gyroradius  and  magnetic  field 

orientation  would  provide  two  additional  parameters,  't'vpical  values  of  M,  kT, 

aiiLi  Ajj  typical  of  the  low  liarth  orbit  environment  are  shown  as  functions  of  altitude 

in  Tabic  1.  '.'alues  of  tile  ambient  electron  Mach  number  M  ,  a  la-  also  shown. 

e 

Table  1.  Typical  Values  of  Linvironmental  Tarameters  as 
Tunctions  of  Altitude,  b  is  altitude  in  kilometers,  Mj,  is 
the  ambient  electron  Mach  number,  M^  i.s  the  ambient  ion 
Mach  number,  kT  is  the  temperature  in  electron  Volts, 
and  A[j  IS  the  Debye  length  in  centimeters 


h  ( km) 

M 

c 

M- 

1 

kT(e\  ) 

Aj^lcm) 

100 

0.  1 

20 

0.  02 

1.  0 

3  00 

0.  Otj 

to 

0.  OMtl 

0,  1  -  0.  7 

aOO 

0.  03 

3 

0.  1:')") 

Ii.  3  -  0.  h 

1000 

0.  02 

4 

0.  25;. 

1 .  0 

The  disk  model  used  in  the  calculations  is  shown  in  1  igurc  1.  .Vole  that  it  is 
not  circula;-,  but  octagonal  and  designed  to  fit  on  a  stjuare  nu-sh.  \Vc  have  chosen 
the  following  sets  of  the  dimensionless  parameters:  ill  \  •  1000,  p  -  lUO,  and 
M  -  -4  and  ii  and  (2)  \  =  200,  p  ■  .aO,  and  M  -  a  (sew  Section  1).  The  \:ilucs  of  \  is 
consistent  '.'.  llh  moderate  eharging,  on  the  order  of  20  to  !  00  \  ,  TvpU'al  \  iJues 
of  in  low  Ikirtli  orbit  are  on  the  order  of  eCrUtniete rs ,  so  that  the  disc  r.adius 
implied  b\  llq.  i2i  is  on  the  order  of  nn-ters.  Mach  -1  .md  !i  ai'e  sonie'Wli  a  less 
til. in  typical  eif  low  p.irtti  oi'bit  s.ilcllites  (on  the  order  of  M.o  h  lOl.  Helow,  we 
diSi  uss  tile  piis.sii.i!  fteitures  assoiiited  with  the  imniericid  results. 

File  qualitalr.e  featiiia-s  of  tlie  numerical  results  can  be  understood  h;isi.(l  on 
.1  u.iive  dis.  iis.sion  of  the  inlt.  r  i e t ion.  The  euises  presented  here  arc  lor  lilgii  voltagc 
ra'g.'.tr.e  ■  Oirging,  ivpic  d  of  .auirging  during  ion  beam  operations.  .As  result  ol 
tae  liign  voltage,  .vtm  (i  l.s  la.'-ge  in  com[)arison  with  the  energies  of  tlu'  ambient  ions 
ind  electiams,  tin-  interiiction  neiir  the  disk  is  qualil:itively  siniil.ar  to  lliat  of  ,i  high 
voltage  probe  in  i  sl.itr  pl  .sni  e  while  the  .iiiibient  ions  .strean.  from  i  iin.  the. 


t  v'"'  ‘ 


'va  X ; 


are  acce  le  ra  teJ  lov.ai\i  llu*  dink  l'i*urr.  all  diri'eliuirs  v.itii  ions  inip.irtui^  tin-  disk  l'r<jn 
holii  ram  ^^nd  AatM.-.  Tiu-  ram  poleiili.»in  ..iid  ipari.hrr  di-nnition  are  >  lose  li- 
exported  in  a  statn  planm.a,  a  ilh  nuitabl*-  nralil  Uatioiin  due  to  the  ran,  motion  oi 
the  ions  >1  sheath  i^>  fornu  vl  and  tin*  potential  'ollov.n  «  -.Mdl  underslooil  (jrolih 
I  See  t'oilov.ing  pa  I’ap- I'aphs ) .  rtif  namber'  ih-nsit’.  d«-..  n-ane.-i  ie .on<U(.>ni i  a !  1  s  I'rorr  lia- 
siieatii  to  tne  vhnk  in  the  n.anru-r  rotiuiia-d  to  .  on~>o:  vc  n,-.-,-..  la-eaunc  oi  tpi 
.‘.rapping  t-h'  tne  ions  ari,,and  ine  dmu,  ihm  re  h*  1 1\ « •  I  •.  uivh  d-uir-.n;.  '  on  iia-  i-rdcr  al 
but  U‘SS  thar'.  ari'.hient'  r«-t.’i‘ai  fAiiM-ioE-  !«■  tne  '-ionnif.  :i-^-iot;  nuiaouiain  tiu-  i.ii.-.k 

in  both  ran.  and  .‘.akt-. 


1';  ir  ;  1  ;•  A  aj-v«-  1 n;  .ia  at  1  i.i  t  i'.  t  i  •-  ami  :  -  ir :  ;  ;  n  .  ;  ■  n.r 

I  mb  lent  pi  I  .-,n:  ;  b  .  t:p  •  di  .ni  •  .a  i  a  .no n  :  <  r  Man  m  :  i  c  . :  <  o  .a  i  . . 
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wake  model  of  Alpert  et  al  also  for  Mach  8.  This  figure  is  not  at  all  fundamental 
and  has  the  status  of  a  cartoon.  The  effect  of  the  geometric  sliadowing  is  to  shield 
the  ambient  plasma  behind  the  disk,  with  the  density  being  heavily  depleted  near 
the  disk  and  approaching  ambient  in  the  far  wake,  t'ombined  with  the  densities  in 
the  near  wake,  the  density  profile  is  expected  to  behave  qualitati'.  els  .is  shown  m 
Figure  3;  a  depletion  required  by  mass  conservation  ne.ir  the  disk,  .1  |)eak  or  knee 
near  one  (ram)  sheath  thickness  in  wake,  a  deep  depletion  due  to  geon.etric  sh.idow 
ing  outside  the  near  svake  region,  and  gradual  approach  to  ambient  detisits  in  tue 
far  wake,  f'igure  2c  shows  the  effect  of  increasing  the  shadow  ing  r. alius  to  tin- 
sheath  radius,  a  concept  to  be  discussed  in  more  detail  in  the  following  p.ir.ig  r.ipii.a 


Figure  2a.  Geometrical  Shadowing  by  a  Disk  .it  Mach  B 


I  igu.'e  2b.  1 'henomenologieal  Inter  action  and  W.ike  Model  to  r  1  lig  a  i 

it  M  o  ll  a.  The  ne.ir  .vake  is  a  i.ae -dimensional  nindi-l  dislorti  ii  to  :  ;  irMi.u.i  r,, 
!i^r..  I'he  wake  model  IS  geometric  shadowing  bv  the  disk.  riu-  i  oUt<  . ;  ^  n  [  ts  ! 
lor.  iiiiii. her  leiisitv;  the  levels  .ire  sep.irated  hv  one-si.xteeutii  ol  tin  .a  hr  a'  a  i’ 


■i .  A  Ipe  j-t .  i  .1 .  I ,  (ill  I'e',  it-h,  \  ,  .mil  1  'itai-VTki  1,  1  ,  1  .  '1  I 
'I'.  itti  A  rtifici.il  Satelliti's,  C 'ons ult.ints  Hure.iu.  New  ii.ia,. 


Figure  3.  Qualitative  Model  of  the  Ion  Density  Profile  Along  the  Ham /wake  Axis 

The  foregoing  discussion  leads  to  the  conclusion  that  peaks  in  the  ion  density 
are  likely  to  appear  in  the  wase,  a  feature  which  is  apparent  in  the  solutions  given 
in  the  following  paragraphs.  .A  simple  argument  indicates  that,  if  they  occur  at  all, 
such  peaks  are  moie  likely  for  high  potentials  and  large  Mach  numbers  than  other¬ 
wise.  The  effect  of  the  high  potential  is  to  accelerate  the  ion  distribution  function; 
with  increasing  velocity  the  ion  distribution  function  narrows  rapidly  (Ciurevich  and 
Meshcherkin^).  The  acceleration  of  the  ramming  ions  into  the  wake  can  thus  be 
expected  to  lead  to  narrow  ion  structures  on  the  order  of  a  sheath  thickness  in  wake. 
Low  potentials  and  low  incident  Mach  numbers  have  the  effect  of  broadening  these 
structures  and  for-  sufficiently  low  potentials  .md  Mach  numbers  they  certainly  will 
not  occur. 

5.  Gurevich,  A.  \' .  ,  and  Meshcherkin,  .A.  i’.  (i!'81>  ion  acceleration  in  an  expanding 
plasma.  Soy.  Fhvs.  .ILTl*,  .“idlNo.  aLhdT. 


For  purposes  of  comparison,  wo  show  in  Figure  4a  a  POL.Alt  contour  plot  of 
the  electrostatic  potential.  This  figure  clearly  shows  the  division  of  the  potential 
structure  into  near  wake  and  far  wake.  To  a  first  approximation,  the  near  wake 
is  symmetric  in  ram  and  wake,  much  like  the  potential  about  a  disk  in  a  non- 
flow'ing  plasma.  The  far  wake  is  a  quasineutral  wake  similar  to  the  geometric 
shadowing  of  the  downstream  flow  by  the  disk.  The  size  of  the  near  wake  in  com¬ 
parison  with  the  far  wake  is  an  artifact  of  the  Fdl.Ali  code.  The  shadowing  of  the 
ions  in  POL.^  H  is  due  to  the  spacecraft  only  and  not  the  spacecraft  plus  sheatli; 
this  problem  is  being  corrected.  Intuitively,  one  would  expect  that  the  narrowest 
part  of  the  far  wake  to  be  about  as  wide  as  the  near  wake.  Figure  4b  shows  ion 
trajectories  from  PC^L.AH;  it  illustrates  the  natui'e  of  the  ion  flow  on  both  sides 
of  th  ■  disk,  as  discussed  above. 

Potential  contours  from  MAt'li  are  shown  in  I  igui'e  4c.  In  contrast  to  the 
POLAR  results,  the  shadowing  in  the  wake  is  determined  by  the  sheath  radius 
rather  than  the  disk  I’adius.  The  line  labeled  "'i'"  Rh.  potential  \  -  1;  the 
remaining  contours  represent  x=  0.  v.a,  0.  li,  ...  64,  LIB.  The  ion  trajectories, 
w'hich  determine  the  charge  density,  have  been  computed  for  potentials  above 
X.  -  0.  25;  at  lower  potentials,  a  geometric.il  shadowing  model  has  been  used.  Tlie 
calculation  was  run  at  rather  tow  spatial  ix-solution,  as  is  evident  from  the  lengths 
of  the  straight  line  segments  used  to  represent  the  contours.  The  accuracy  of  tlie 
contours  in  the  tar  wake  is  not  known  at  the  present  time;  it  is  possible  that  thev 
will  smooth  out  with  further  relaxation  of  tlie  solution.  The  fact  that  the  edges  of 
the  far  wake  as  determined  by  the  \  ^  0.  25  contour  are  parallel  is  an  artifact  of  the 
method:  the  boundaries  ol  the  computational  mesh  jirevented  it  from  exptmding 
further.  Llectron  number  density  contours  are  shown  for  the  same  solution  in 
Figure  4d. 

:i.  COMP AKl.SON  Ol  MACH  W  PI  II  SIMPl.i;  ANAI.YTIC.AI.  MOOCL.S 

The  purpose  of  this  section  is  to  compare  M.ALll  with  simple  an. dvtical  models. 
The  models  used  in  ram  and  wake  are  reasonably  good,  while  the  model  used  .ilong 
the  edge  of  the  disk  is  entirely  qualitative.  The  purpose  in  these  comparisons  is 
to  est:iblish  a  i-e.isonable,  if  inexact,  basis  for  judging  the  accurac>  of  the  moilels. 
The  physical  mech.inism.s  invoUed  m  the  spacecraft  interaction  are  somewliat 
more  complicated  than  the  simple  models.  The  extent  of  the  agreement  outlined 
in  the  following  [i.iges  suggests  th.»t  liie  models  jierform  reasonahlv  well. 


MACH;  ELECTRON  DENSITY  CONTOURS 


\f/  =  -2e4-/(m\'  )  =  M" 


4  =  (p/M)x/a 


and  is  the  location  of  the  sheath  edge,  while  the  equation  for  tile  numbei-  density 


n  /  n  =  1  /  N  (  1  -f  v ' 
o 


where  is  tlie  ambient  ion  number  density.  In  figure  5,  tlie  data  points  an-  the 
values  computed  by  MACH  while  the  curves  are  computed  from  liq.  (4>.  The 
agreement  of  the  values  is  everywhere  within  a  few  percent  of  the  pe^k  potential. 


HACH  il:  ^  POTENTIAL.  A  DENSITY 
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figure  5.  Comparison  of  MiXCIl  Potential  and  Ion 
Number  Densities  in  Kam  for  Mach  4  and  8  Witli  Oiu 
Dimensional  Sheath  Contours 
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U  should  be  emphasized  that  there  are  no  xi-ee  parameters  in  tlie  analytic  solutions; 
for  the  given  plasma  conditions,  the  solutions  are  determined  completely  by  the 
potential  on  the  object  and  by  the  requirements  that  the  potential  and  eli  ctru  field 
both  vanish  at  the  sheath  edge. 

Figure  t>  shows  the  comparison  o.  M.XC'll  potential  and  number  density  profiles 
.dong  the  rim  direction  for  Mach  8  with  an  analytic.il  approximation.  I’he  rim 
direction  is  here  defined  to  extend  from  the  disk  along  ..  line  pe rp  'iidicnJa r  to  the 
ram-wake  direction.  The  approximation  is  not  vers  good  .md  is  given  for  purpose.-, 
of  qualitative  comparison  only.  In  tliis  approximation,  the  potenii.il  is  computed 
assuming  that  it  is  cylindrically  symmetric  around  the  edge  of  the  disk  .ind  th.d  the 
ion  number  flux  into  the  symmetric  sheath  is  determined  by  the  ion  therm.d  speed 
and  not  its  bulk  speed.  The  curves  are  the  solutions  of  I'oisson's  equation 


(l/^)d(4dx/d4)/d^  =  2^  /l^\ll-x\il 


I  81 


n/riQ  =  n+  X)] 


I  h) 


I  =  ‘'/A, 


I  I  01 


along  with  a  phenomenological  modification  |r  -*  \(r“4  d"i,  uhere  d  0.  2a  is  the 
thickness  of  the  disk]  designed  to  represent  the  potenti.d  .is  it  approaches  .i  disk  of 
finite  thickness.  It  is  seen  that  the  potentials  are  in  reasonably  good  agreenietii. 
within  a  few  percent  of  the  peak  potential.  .Again,  it  is  to  be  emphasi.’ed  th.d,  .iside 
from  the  construction  of  the  model  (which  is  only  a  crude  represent.dion,  at  hesti, 
there  are  no  free  parameters. 

Figure  7  shows  the  comparison  of  MACH  potential  and  ion  number  density  pro¬ 


files  on  the  wake  axis  for  the  same  Mach  numbers  with  .in  analytic  .solution  gum 
by  .Alpert  et  al;  ^ 


2  2  2 

n  =  n  exp(-M  a  /z  l 
o 


1  1  t 


■.'.here  z  is  the  distance  from  the  disk.  The  agreenu'iit  between  cum[iutation  .md  Uie 
geom.etric  shadowing  model  shown  .is  solid  lines  i.s  re.tsonablv  gwod  m  the  t  .r  '.'..me 
region,  beyond  perhaps  thiee  to  six  disk  raoii.  Phe  overshoot  ul  eight  disk  radii 
for  Mach  4  and  t'.velve  disk  radii  for  M.icli  H  is  likelv  ol  numern.d  origin,  il  .so,  it 
should  be  possible  to  cure  it  by  changing  the  mesh  size.  It  should  be  noted  th.d  ttiis 
type  of  agreement  is  not  expected  to  hold  gcnerallv,  but  is  a  i  esult  of  tlie  f.n  t  tli.d 
the  case  presented  can  be  adequately  represented  by  a  quusineutr-al  far  iv;ike.  If 
we  h.id  assumed  that  the  electrons  were  dense  and  energetic,  iiith  perhaps  the 
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HACH  8:  POTENTIAL. 


DENSITY 


DISTANCE 


I'lfTure  6.  Comparison  of  MACH  Potential  and  Ion  Number  Densities  Along  the 
Itim  Uireetion  \Vith  Cylindrical  Sheath  Contours 


10  kc\  energies  characteristic  of  the  auroral  environment,  we  would  expect  tfiat 
the  structure  of  the  far  wake  would  be  largely  controlled  by  ambipolar  diffusion- 
t!ie  wake  should  he  substantially  shortei'  because  of  tlie  transverse  acceleration  of 
ions  by  thi-  diffusion  process. 

The  daslii’d  lines  in  1-  igure  7  show  the  effects  of  increasing  the  I'ffective  shadow¬ 
ing  radius  by  25  [lercent  foi-  Mach  -1  and  10  percent  for  Mach  8.  The  physical  bases 
IS  that  at  least  [)a  i  t  of  the  sheath  must  be  responsible  for  geometric  sliadowing. 

These  increasi-s  constitute  about  one-lnlf  and  one-quarter  of  the  rim  sheath  radius, 
re-spectively .  \‘vhile  increasing  the  shadowing  radii  to  tlic  full  slieath  radii  gives 
very  bad  fits,  the  smaller  increases  can  reasonably  account  for  the  far  wake  de-isity 
within  a  few  pi-rcent  of  .mibient  beyond  about  three  disc  radii,  with  tin-  exceptions 
of  the  overshoot  r-egions.  This  gooil  agi-eement  suggi-sts  that  tlie  i-ffei  t  of  the  sheath 
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The  discussion  so  I'ai-  lias  centered  on  tlie  analytic  comparisons  ishich  it  is 
possible  to  make.  The  M.TCH  computer  model  is  in  good  agreement  with  tlie 
anals  tic  tests  witli  the  exception  of  the  ion  density  in  rim.  .\side  fi-om  the  rim 
densitv,  the  outstanding  feature  which  does  not  fit  an  analytical  model  i.s  the  e.xi.st- 
ence  of  the  peak  in  the  ion  number  density  inside  about  three  disk  i-adii  in  uake  lur 
Mach  h.  Wliile  we  have  not  established  that  it  is  a  genuine  feature,  e  believe  this 
pe.ik  IS  likely  to  be  real  based  on  the  discussion  associated  with  l  igiii  es  2  and 

3.  The  peak  follows  the  general  qualitative  trend  outlined  e.i  idler,  UK  luding 
appearance  close  to  the  disk  and  at  high  Mach  numbers. 

4.  COMH.MUSON  OK  M.ACIl  .AM)  I'OL.VR 

To  test  the  1HJI..-\.H  code  against  MACH,  we  have  run  I’tlL.AH  for  the  same 
conditions  mentioned  .above.  In  this  section  we  compare  HOI..\H  results  with  tiiose 
of  M.AC'II.  The  bidef  discussion,  which  covers  some  of  the  material  discussed 
above,  is  intended  to  address  the  relative  accuracy  of  both  models,  f.stiniaies  of 
accuracy  are  in  no  way  absolute  but  represent  a  judgment  based  on  comp.irison  of 
the  models  with  each  other  and  with  analytic  approximations.  It  should  be  noted 
'‘  at  all  of  these  comparisons  are  at  Mach  5  rather  than  at  Mach  4  .ind  li. 

f  igures  tia  and  lib  show  POL.AK  and  M.-Xt'll  profiles  of  electrostatic  potential 
and  ion  number  density,  respectively,  in  ram  for  Mach  5  along  with  the  solutions  of 
Llq.  (41.  The  i’OL.XH  potential  is  too  broad  and  demonstrates  a  known  tendency  of 
POL. .AH  to  increase  the  slieath  thickness  by  about  one  grid  spacing.  .As  mentioned 
earlier,  this  problem  can  presumably  bo  controlled  by  r-unning  at  higher  I'esidution: 
the  operational  problem  is  to  balance  increased  aiwuracy  against  incre.ised  eompii- 
tation  time.  The  POLA  H  densities  are  essentially  exact  in  comparison  with  l  .q,  (  41 
till'  .M.AC'll  densities  are  somewhat  too  high,  a  feature  not  seen  in  the  Mach  1  and  K 
results.  These  ei-rors  suggest  that  the  errors  in  i'f)l,.AK  and  ,M.'\t'll  a  I'e  cui-rentlv 
on  tl'.e  order  of  10  percent  of  peak  values  in  ram  at  the  spatial  resolution  shown. 

figures  Ua  and  ''b  .-.how  the  same  profiles  in  the  rim  direction.  rtie  s.une 
teiidencv  toward  ve.  i,- la-xpansion  of  the  sheath  is  .ipparent.  The  densitv  profilca 
show  f luctu.itinn.s  on  the  order  of  10  to  20  percent  of  ambient  which  mav  be  i.aken 
III  be  t.'ie  order  of  ni.ignitude  of  the  error  in  the  models  at  low  resolution. 

1  igu.i’e  10  -.iio.i.s  lull  numbe  r  dens  1  tv  contours  in  wtike.  tleometi'ii  ah.ido  a.  ing 
b.  tne  disk  fo;-  M  icli  "i  is  shown  .is  a  dashed  line.  M.ACll  shows  the  geonieti  i, 
anolowing  bv  tne  disk  mentioned  .hove.  riu-  ditleri'nee  iH.'tween  M.AC'll  .md  geo¬ 
metric  shadowing  .ippears  to  be  unphysie  il:  the  M.Ac'H  data  correspond  to  geonietri 
sh.idowing  hv  .i  disK  smaller  than  the  actual  disk  si/.e.  P(l|..AK  shows  similar 
sliadowing:  the  displacemi.-nt  could  reasonablv  be  interpreted  ..s  shadowing  due  to 


disk  plus  sheath  except  lor  the  fact  that  only  pure  disk  shadowing  exists  in  I'Ol.Ait. 
As  such  attribute  errors  on  the  order  of  10  percent  to  20  pei-cent  to  both  WAt  11 
and  POLAR  beyond  a  feu'  disk  radii  in  wake.  Both  codes  show  lai-ge  lensitv  [H-aks 
in  wake,  M.\C11  to  twice  ambient  and  POl-.vK  to  about  0.7  time.s  anibu  nt.  Both 
figures  are  surprisingly  lugh.  The  MAC.'H  density  peak  appears  to  be  unreali.-,ticalh 
high.  While  we  have  argued  that  density  peaks  in  wake  are  likelv  to  be  real,  ihe.''e 
large  enhancements  leave  open  t.he  possibility  that  errors  in  the  near  wake  .ire  on 
the  order  of  100  percent,  l  >nly  lurther  compai'ison  ..nd  jjossible  rraniifu  ation  of 
botii  codes  will  eslablis'i  the  ii-ulh. 

->.  DISCLSSION 


1.^  e  h.e.e  coni[>e.red  the  results  of  the  M.Ai'II  and  l’(ll..\R  computer  codes  fo  r 
conditions  Ispic-L  of  n.oderate  lu  high  voit.ige  neg.or.e  charging  for  satellitc.s  m 
low  Lartii  oi-hit.  We  find  that  the  m.odels  a.pii.',.:-  !o  be  in  genei'.il  quantitati'.  !■  .iiiu 
qualitative  agreen.ent  with  Ih.e  features  of  the  interaction  whlcil  are  susceptible  to 
.jheck.  Tlie  ma:or  exception  is  die  tendency  of  tlie  POL.AR  wake  to  la  p  i  losent 
shadow  ing  bs  t.he  disk  rather  tlian  the  sheath.  V.  e  reiterate  lira  tile  .nialvtic  ciiecks 
in  ram  and  r'lm  have  nf>  free  [la  ramete i-s  .aside  f.-'orn  the  si.-ui  turt  of  the  models. 
The  agreement  between  the  i-odes  and  tne  .inalvtu-.d  n.udel.s  in  ihese  regions  range 
from  reasonablv  to  I'emark-ibly  good,  V\  hile  the  agreement  in  the  f.,  r  wake  has 
one  free  parameter  i  the  efiective  geometric  shielding  radiusL  the  n.iture  of  the 
shadowing  is  in  i-easonable  agreement  with  the  phy.sics  of  the  interaction.  In  addi¬ 
tion,  We  have  argued  that  the  occurrence  of  ion  density  iietiks  near  the  disk  in  ,ike 
IS  physically  reasonable,  at  least  for  high  eiKJugh  potenti.i!  and  Mach  rnimber. 

Based  on  this  analysis,  we  believe  that  the  MAt'l!  and  I'OL.AR  computer  codes 
adequately  compute  the  spacecraft  interaction  and  wake  structure.  J-'urther  work 
is  required  to  resolve  the  disagreement  between  PtJI.AH  and  M.Af  ll  on  the  peak  ion 
density  near  the  disk  in  the  wake. 
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